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Abstract 

An ASEP with two species of particles and different hopping rates is considered on a ring. Its 
integrability is proved and the Nested Algebraic Bethe Ansatz is used to derive the Bethe Equations 
for states with arbitrary numbers of particles of each type, generalizing the results of Derrida and 
Evans [10] . We present also formulas for the total velocity of particles of a given type and their limit 
for large size of the system and finite densities of the particles. 
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1 Introduction 



An idea which has proved to be quite useful in understanding the behavior of systems out of equilibrium 
is to study solvable models. An example of such models, which has driven a lot of attention for at last 
two decades, is the asymmetric simple exclusion process (ASEP) Tj. It describes a driven lattice gas 
[2 [3] where particles can hop on adjacent sites with asymmetric rates and hard core exclusion. 

Different methods have been applied to study the ASEP, and each of them seems better suited to 
study certain aspects of the problem. The Matrix Product Ansatz for example has been employed with 
success in determining the density profile of steady states, steady currents or diffusion coefficients (for 
reviews see [4J [5j E] ) . On the other hand other quantities like the relaxation time are more easily dealt 
with by means of the Bethe Ansatz [7J [5] . 

Actually in [S], Derrida and Lebowitz showed how a modification of the Bethe Ansatz of Gwa and 
Sphon [8 , could be used to compute the full large deviation function of the time averaged current for the 
ASEP with one specie of particles. 

Shortly later Derrida and Evans [lOj considered the problem with a second specie of particles, and 
thanks to a Bethe Ansatz they were able, not only to reproduce known results [TT1 [T2"] about the phase 
diagram of the steady current of a particle of second type, but also to compute its diffusion coefficient 
and in principles all the higher cumulants. 

Multi-species generalizations of the ASEP have been considered in several papers [13l HH [HI HS1 HZ] ■ 
The fact that they are integrable is not at all obvious. The most natural integrable generalizations of 
the single specie have a hierarchical structure based on quotients of the Hecke algebra, which ensure the 
integrability, as explained in [T8] . 

On the other hand, in the ASEP with two kind of particles and different rates, as considered by 
Derrida and Evans, the hierarchy is partially spoiled precisely by the different hopping rates, and to 
understand its integrability from a point of view of the Yang-Baxter equation one cannot make resort to 
the Hecke algebra commutation relations. 

Our first point in the present paper is to make manifest the integrability of the ASEP with two species, 
by showing an R— matrix which solve the Yang-Baxter equation and gives, through the usual procedure, 
the transition matrix of our problem. Once we have this we employ the machinery of the Algebraic Bethe 
Ansatz (ABA) to derive the Bethe equations and the eigenvalues of the transition matrix (see [T!5] for a 
review on ABA, and [20] for recent application of ABA to the ASEP). 

Since we have a number of species greater than one we are led to perform a Nested Bethe Ansatz. In 
the case of arbitrary number of species, but with hopping rates independent of the types of particles the 
Nested Bethe Equations have already been derived in [2T] . 

With the Bethe Equations at disposal we can tackle the problem of determining the cumulants of the 
total velocity of particles of a given type, or of joint cumulants, in presence of an arbitrary number of 
particles of each kind. We present the exact formula for the average velocity of particles of second type 
and consider as well the limit of large size of the system with finite non-zero densities of particles. We 
comment also on the difficulties about the determination of the higher cumulants. 

The plan of the paper is the following. In section 2 we show the integrability of the ASEP with generic 
rates by presenting an R— matrix which solves the Yang-Baxter equation and generates the transition 
matrix of the ASEP. In the same section we use the techniques of the Algebraic Bethe Ansatz to diag- 
onalize the transition matrix, arriving at a set of Nested Bethe equations. In section 3 we analyze the 
Bethe equations and derive the exact formula for the total velocity of M-i particles of type 2 on a ring of 
size N and in presence of M\ particles of kind 1, we comment on the derivation of the higher cumulants. 
The large N limit of the velocity is worked out in section 4 where we show that for non-zero densities 
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of particles of each type, there is no phase transition. In Appendix A we sketch the derivation of the 
Bethe Equations for an ASEP with twisted boundary conditions, that we need for the derivation of the 
Nested Equations. In Appendix B we study the integrability of models with higher number of particles 
and arbitrary hopping rates. 



2 Yang-Baxter for two species and different rates 



In [10], Derrida and Evans have employed the coordinate Bethe Ansatz to study an ASEP in presence 
of an impurity, which in the following will be treated as a second specie of particles. If we indicate with 
the empty site (which can be considered as a particle of type zero), 1 a particle of first kind and 2 
a particle of second kind, then the rules that govern the stochastic evolution of the system during an 
interval of time dt are purely local on couples of neighboring sites and are given by 



10 
20 
12 - 



► 01 with rate 1 
• 02 with rate a 
21 with rate f3. 



The fact that the problem can be solved by Bethe Ansatz, as done in [10 for the case of a single 
second type particle, means that it is integrable. Our first task is to understand better its integrability 
showing the Yang-Baxter equation behind it. 

The transition matrix of our system can be written in terms of matrices that encode the local transition 
of particles. Let us define the basis of the local space of states as: 

|0) = empty = particle of kind 0, 
|1) = particle of kind 1, 
|2) = particle of kind 2. 

In the basis (|0), |1), |2)) a ® (|0), |1), |2))f, we introduce the matrices 
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The ASEP with two species is defined by the following equation for the probability of a given configuration 
C 

j t P t (C) = J2 M °' ' (C,C')Pt(C'). (2) 
c 

where 

M UW ' ,/2 °' V12 {C,C') = ^(^ (10) + uEf 0) + (3El 12) ). (3) 

i 

Since we are in presence of several species of particles we can introduce the relative distances covered 
by particles after an initial time t = 0. By this we mean the distance Y l ° covered by all the particles of 
kind i with respect to particles of kind j. It increases of a unity each time a particle of kind i jumps to 
the right of a particle of kind j, and decreases of a unity when the opposite happens. In our case we have 
three kind of particles 0, 1, 2, hence we consider Y w , Y 20 , Y 12 . The joint probability P t (C, Y 10 , Y 20 , Y 12 ) 
of being in a configuration C, and having Y t y — Y y satisfies an evolution equation which is better written 
in terms of a generating function 

pv 10 ,V2 ,vi2(Q^ _ e "ioY 10 +v 20 Y 20 +v 12 Y 12 p^Q ylO y20 y 12 ) (4) 

ylO y20 yl2 

The evolution equation satisfied by F" w ' V20 ' Vl2 {C) is 

^_p v Wy via,v\i _ ^2 M Ul0 ' l ' 20 ' u ' 12 (C Q'~ s jP Ula - U2a - Ul2 (C') (5) 

dt c , 

One obtains {e VwY ? '+ V ™ Y ? '+ v ^ Y t'' ! ) summing F" 10 '" 20 '" 12 (C) over C , hence its large time behavior is 
determined by the largest eigenvalue \(vio, V20, v \i) of the transition matrix M Vl0 ' V20 ' Vl2 (C, C) 

^wY t 10 +u 20 Y 20 +u 12 Y t 12 ^ ^ e X{v 10 ,u 20 , Vl2 )t_ 

We find such an eigenvalue by employing the Algebraic Bethe Ansatz. 

Our first step is to find an R— matrix which satisfies the Yang-Baxter equation, the inversion relation 
and such that its derivative reduces to the linear combination of E^'s matrices 

E^ + aE<M + I3E {12) 

Once we have this we construct the transfer matrix in the usual way as trace of product of L— matrices 
(L = PR) and we are insured that its logarithmic derivative will be the desired ^ i (i?j 10 ^+a£ , l - 20 ^+/?£'j 12 ^) 
We provide a solution of the Yang-Baxter equation 

R a ,b{y, z)R b , c (x, z)R atb (x, y) = R b , c {x, y)R a ,b{ x , z )RbAv-, z )- ( 6 ) 

of the form 

R(x,y) = l + g w (x,y)E (10) + g 2 ^y)E^ + g 12 {x,y)E^ 2 \ (7) 

where 

9i2(x,y) = l- ^—±——± -, g w ( x ,y) = l-e* V; g 20 (x, y) = 1 - ^—-i-—^. (8) 
We define the monodromy matrix of a system of size N as 

^a®Jif(x,f)) = L a . aN (x,r] N ) . . . L ata2 (x,r]2)L atai (x,r]i). (9) 
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where L a ^{x,y) = P a .bRa,b{x,y), and P a ,b is the permutation operator, i.e. Pv a (£> Vb = Vb ® v a . . The 
transfer matrix is given by 

T(x,ff) = tr a ^a®Jt?(a, ff) 
Thanks to the Yang-Baxter equation ^ we get 



^a®j#>(z, y)-%®j?(x, y)R a , b {x, z) = R a ,b{x, z),9'a®^{x, y)%®j#>(z,y). 



(10) 



and, tracing over the auxiliary space, we obtain that the transfer matrices with different values of the 
spectral parameter commute among themselves 



[T(x,ff),T(x',ff)]=0. 



(11) 



The transition matrix of our system is obtained choosing r\i = 0, and taking the logarithmic derivative 
of T(x, ff) at x = 



M 



^10^20^12 



(C,C) = -T(0,0) 



_idT(x,Q) 



dx 



x=Q 



(12) 



We come now to the Yang-Baxter algebra, which can be easily read from eq. (|10|) . Let us write the 
monodromy matrix as 

A(x) Bi(x) B 2 {x) 
■K = | Cx{x) D n (x) D 12 (x) j (13) 

C 2 {x) D 21 (x) D 22 (x) 



The transfer matrix can then be written as 

T(x) = A(x) + D n {x) + D 22 {x), 



(14) 



where we have simplified the notation omitting £ which is fixed to be zero. For later purposes let us 
rewrite the R— matrix as 



R{x,y) 



( 1 

1 (1 

1 
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e x-V 














R (l) 

-"-11,11 

(1) 

12,11 



(1) 

21,11 



(1) 








p(l) 

-"-11,12 

(1) 

12,12 



(1) 

21,12 

(1) 



(1- 






l+q(e*-lU „ 
l+o(e»-l) I 





l+q(e x -l) 
l+a(eB-l) 





\ 



R 



-"11,21 

(1) 

12,21 



p(l) 

-"21,21 -"21,22 



p(l) 
-"11,22 

-"12.22 



•■22,11 -"-22,12 

The matrix R^\ which actually depends on x and y, is given by 



RW(x,y) 



( 1 



(1 

V o 





l+g(e-»-l) 
l+/3(e— -1) 
l +| a( e -«_l) - 

l+/3(e- a: -l) ^ 







\ 



1 
1 / 



R 



(i) 



R 



(i) 



(15) 



(16) 



It is nothing else than the R— matrix corresponding to the ASEP with a single specie presented in eq. ([63 
in appendix A, with a different parameterization 

1 



Pie- 



ty 



(17) 
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With this notation we can write the commutation rules of the operators A, Bi, Ci and Dij appearing 
in i^, for different values of the spectral parameters 

[A(x),A(y)]=0; (18) 

A(x)B 1 (y) = e „ 10 (f, _ eV j B\ (y)A(x) + ^^-^ ) B 1 (x)A(y); (19) 

A(x)B 2 (y) = e lt %Z% B 2 (y)A{x) + J+^Z^ B^Aiy)- (20) 

Bi(x)Bj(y) = B l {y)B k {x)R^ lk {x,y); (21) 

Dy(x)fl fc (i,) = e „ 10 g_ ea) (B m (y)D lB (z)i?g m „ (a, y) - B j (x)D lk (y)); (22) 

D 2j (x)B k (y) = e ffigl£) (B m (y)D 2 n(x)R^ jk ,mn(x, y) - B J (x)D 2k (y)). (23) 

The ansatz for an eigenvector of the transfer matrix, keeping into account the non commutativity of 
the BiS for different i, is given by 

\* Mi ' M °(yu...,vr))= E Kl:t B M--- B M\\i), (24) 

ii,...,i r 

where ||1) is the reference state, defined by 

/ 1 






V o 



and it is an eigenstate separately of A(x), Dn(x) and D 22 {x) 

A{x)\\l) = \\l), D ll {x)\\l)=e N ^{l-e x ) N \\l), 

D 22 (x)\\l) = (ae^) N (l-e*) N \\l), 
and correspond to a completely empty system. The labels M\,M 2 in eq. ((2^|) mean that we require M\ 



Bs of type 1 and M 2 Bs of type 2, i.e. we are restricting to the sector with Mi particles of type 1 and 
M 2 particles of type 2. 

The eigenvector equation for \^ Ml ' M2 (yi, y r )) puts constraints on the j/jS. Let us first apply the 
operator A{x) to \^ Ml,Ma (yi, ■ ■ ■ , y r ))- We get a wanted term, i.e. a term proportional to the vector we 
start from, of the form 

e \ I 1 + a(e x — 1) \ -p-r 1 



and unwanted terms of the form 



II j-^l^Hyu...^)) (25) 



e*\ Ml /l + a(eW-l)^ M2Ml+M2 



IT ( pVl _ P y,) B ^ ® ■ ■ ■ ® B iVj-i) 



M( yj ,y)M( yj ^,y) . . . M (y u y)^ M - (y u . . .,y r )\\l) (26) 

with 

M(y 3 ,y) = i?,g 2jiJ ,( % , yj> i)i?g 3j2j , ( yj , %+2 ) . ■•^ lirjip _ lJr (%»yi-i) 

The unwanted terms have to cancel with similar terms coming from the action of Dn(x,y) + D 22 (x,y). 
From the action of D kk we get a wanted term 

Mi+m 2 

w k (g)(l-e x ) N II 7^^B( yi )^-.-^B(y Ml+u jTg\x,ff)^ Ma ( yi ,..., yr )\\l), (27) 

i=l ^ ' 
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with 

M1+M2 , v v Mi+M 2 , , . 

n (£). *»-^-«» n i 1 ^ 1 

t=l v 7 i=l v 

and 

T^(x,ff)=L^ Mi+M2 (x,y Ml+M2 )...L^ i (x,y 1 ) (28) 

is just the monodromy matrix of TASEP with a single species as explained in the appendix [XJ We get 
also an unwanted term 

n _ eVj) N 1 

^ e *- e v S n (e „ f _ eW) fl w ® ® ■ • ■ ® B (gi-i)x 

M(y 3 , y)M(yj-i, if).-. M( yi , y)* M > (y u ..., y r )T$ {y 3 , y)* Ml < Ms (y u . . . , y r ) | 

In order to get the cancellation of the unwanted terms we first have to diagonalize u/i(y)T^ + ^(ifjT^ ■ 
This is the transfer matrix of a TASEP with a single specie and twisted boundary condition and can be 
diagonalized by the Algebraic Bethe Ansatz as in the case of non-twisted boundary conditions, as done 
in |20j . we briefly recall how it works in appendix lAl Here we simply apply the results explained there. 
One has only to be careful in translating the parameters e Vi appearing in the appendix, following eq. (|17| 

«" 1 (29) 



l + /3(e-»-l) (1-/3)15-1' 

where we have defined 15 = 1 — e Vi , and remember that we consider only the sector with M2 particles in 
a system of size N = Mi + M2 ■ 

For the auxiliary spectral parameters Zs we get the Bethe equation (|74[) which now reads 

e*. a \ * (l-aY^m-l) *f / Zj \_( e^e-o \ MM 

#») l\ il-Yim-l-ZjfXi-l))^ Z k )-{ a e- ) ■ W 

While the cancellation of the unwanted terms coming from A and D^k leads to a second Bethe equation 
The eigenvalue of the transfer matrix can be read from eqs. (|25l [27]) and eq. ([75| 



\ / \ 7 i=l v y i=l i=l v y 

(32) 

Taking the logarithmic derivative in x = we get the eigenvalue of the transition matrix A 

1 dA(x) M1+M2 



M 1 -aM 2 + i- ( 33 ) 



15 



A(0) dx 

3 Analysis of the Bethe equations 

For convenience of notation we divide the Ys in two sets, Y^ with i = 1, . . .Mi, and Y^ with i = 
1, ...M2. The solution of the Bethe equations wanted behaves in the limit — * as l^ 1 ' — > 1, 
1G — > l/a and — > z£ , where the have to be determined. Actually we will see that the Z^s 



(i 



depend on how the limit is taken. For this reason we will redefine i/y — > vvij and take the limit v — > 
keeping fixed. What happens is that the Z^P's depend on these z/y (or more precisely on their ratios). 



From eqs. p0|3ip we get {& 



Mi M 2 
e V2< ) M 2 +vioM 1 a M 2 TT}^ 1 ) TTY.^ 



1 and by continuity 



(34) 



i=i 



Let us introduce the following auxiliary variables 

Mi + M 2 

Q — e (fi2+^20-^io)M 2 +^ioA r J~J (\ — Yi) K = 

i= 

The Bethe equations become 



e » 12 (Mi+M 2 ) Y[t=i (! - Yi) Uit\ z k 



i=i 



(ae 



u 20 -u 10 \N-Mi-M 2 -r-rMj+M; 



Afi+M 2 
i=l 

M 2 

(1 - y,) Ml (l - aY,-) M2 (Mj - l) Ma = Ca M *Yf J[ (6Yj - 1 - Z k (Yj - 1)) . 

fc=i 



(35) 



(36) 



(37) 



We notice that if we keep K as an unknown, combining eq. (|34p with eqs. (|36| l3T |) we recover the definition 
of K given in eq. ([35| . Hence from now on our basic equations are (|34p. (|36|) and (|37| . Following steps 
similar to the ones in |10j we obtain the following representation for the eigenvalue, which generalizes 
eqs. (33, 34 and 36) of [TU] 



E 



dy_l 
2iri y 



where for us 



Q(y) 



1 Jl/a 

y N a M *ni=i(by-i-z k (y-i)) 

(1 - ay) M * (by -1) M *(1- y) M ^ 



Taking the logarithm of eq. l|3"4"|) we get 

v w Mi + v 20 M 2 



00 /^n 

E — 



1/c 



Zm y 



while eq. (j3l)|) becomes (after having taken the logarithm) 

M 2 log(-Zj) = log if + 2i7rj + Mi log(6 - 1) + M 2 log(6/a - 1 - ^-(1/a - 1))- 



(38) 



(39) 



(40) 



(41) 



E 

71=1 



1/a 



6- Zi 



2ni by — 1 — — 1) 



Taking the logarithm of the first equation in (|35|) we get also the equation 



vio(N - M 2 ) + (z/ 2 o + v 12 )M 2 = - 



C" 



71=1 



1/a 



dy 1 
27ri y — 1 



[Q(y)] n - 



(42) 



Now, if we redefine i/y — > i/J/y, we see that eqs.((3Hl SHI HH and give in an implicit form the power 
expansion in v of A. What one should do in principle is to expand log(K) and Zj in powers of C, use 
eq. (|4*Tj) and a combination of eqs. (|4QI42[) to derive the n-th order term of these expansions in terms of 
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lower orders terms, and then work out the expansion of A in powers of v. This way one would find the 
cumulants of the total number of particles flown viqY^ 1 ^ + V2oY^ 20 ^ + v^Y 1 " 12 ^ 

M^oY^+u 20 Y^+u 12 Y^) ) 



X(u) = lim 

t — >oo 



log(e* 



t 



E 



((^oy( 10 ) + ^oy( 20 ) + ^ 12 y( 12 ))») 
t 



c v n . 



(43) 



Concretely this is of course quite laborious and one doesn't find any illuminating formulas in general, 
but one can easily find at least the velocities. Let us work out explicitly the average of the total velocity 
of the particles of second type. For this we have to chose v\§ = 0, ^20 = —^12 = 1, and the velocity is 
given simply by the linear term of the expansion of A in terms on v. In the limit v — > 0, the ZjS satisfy 
a very simple equation 



b 

1 

a 



zf> \ - 

a 



1 



1 



Mo 



whose solution is 



7(0) 



6) e T&F[X(0)(5_ l)Mi]l/M 2 



(44) 



(45) 



The Z^ are now expressed in terms of a single unknown which is determined taking the first order 
in C of the constraint equation |42|) 



l/a 



dy 



2m(y-l) 



)(0) 



(y) = o, 



where is the value of Q for v = 0, which is given by 



? (0) (y) = 



»M 2 



(l-y) Ml \(l-ay) M 3 



_ K (o) 



0) (b 1) 



M1+M2 



(by - l) M * 



Then for we find 



K (0) 







dy 




y N a M 2 


k/a 


2ixi 


(1 


-y) M i + 1 (l-ay) M 2 






dy 




y N (b-l) M i+ M 2 




■fl/a 


2-Ki 


(1 


-y) M i + 1 (by-l) M 2 



(46) 



(47) 



(48) 



Notice that, as stated before, the value of K^°\ and hence of the Q s, depends on the choice of v^. Had 
we chosen different values for v\q, V20 and i>\2, we would have found a different value for K^°\ 

Now we have all the ingredient we need to derive the velocity of the particles of kind 2. We have 
simply to consider the linear part of eqs. ([38"l l40 f 



v 2 = lim = M 2 



dy 



2-ni y 2 



dy QWfe) ' 
2tti y 



§1 + fl/a 

which is more conveniently written in terms of the following two auxiliary functions 



F 



N,M U M 2 



r N.,M ± ,M2 



+ 

J J l/a 



dy 



2tt (y - l) M i (ay - 1) M * 



l/a 



dy 



2ir (y-l) M ^(l-by) M * 



(notice that the contours of integration are the same for the two integrals) 



v 2 = M 2 



rpa rrib jpoi jpb 

r N-2,M 1 ..M 2 T N,Mi + l,M 2 r N,Mi +l,M 2 r N-2,M±,M 2 



r N-l,Mi,M 2 N,Mi-\-l,M 2 



r N , A h + 1 , M 2 r N - 1 , A i 1 , M 2 



(49) 

(50) 
(51) 

(52) 
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When a — I and 6 — we can compute both F a and F b exactly [f| 

*n7mi ,m, = ( Mi + m 2 - l) ' F W,m 2 = ( M *_ 1 ) ■ t > 



And for the velocity we get 

N — Mi - 2M 2 , , 

Wa - Af 2 ^ (54) 



4 Large N limit of the velocity 

We want now to consider the limit N — > oo with non-zero densities of particles of both species p\ = M\ /N 
and p 2 = M2/N. To find the asymptotic formula for the velocity we have simply to determine the 
asymptotic of F a and F b , which are easily given by the steepest descent method. Both integrals have 
two saddle points which correspond to the solutions of the equation 

1 pi p 2 n 



y y - 1 ny-1 

where k = a for F a , and n = b while for F b . The expression for the saddle points is 



(55) 



± _ K + l- px- np 2 ± y/jn + 1 - p! - Kp 2 ) 2 - 4k (1 - pi - p 2 ) , , 

2k(1 - pi- p 2 ) 

It is easy to realize that both saddle points are on the real line, one is situated between 1 and 1/k, the 
other is situated to the right of 1/k. 

When we compute F a we can merge the contour around 1 with the one around l/a and let the 
resulting contour pass through y+, hence F a is dominated by the contribution from y+. In computing 
F b the integral around l/a gives no contribution because the integrand is holomorphic there. The contour 
integral around 1 can now be deformed only to pass through , because of the singularity present in 
1/6, hence F b is dominated by the contribution from y^ . 

In conclusion we get 

v 2 - M 2 C-L + -L - lV (57) 

\Va y b J 

By the same procedure one can find the total velocity of the particles of kind 1 

vi = -{N- Mi) ( -1 + 4 - lV (58) 

y»y b \yZ y b ) 



From the eq. ([56|) we see that the non-analyticities of the total velocities (|57l I58|) are located at the zeros 
of the square root, i.e. at 

(a+l-pi- ap 2 f - 4a(l - px - p 2 ) = (59) 

and at 

(b+l-pi- bp 2 ) 2 - 46(1 - pi - P2 ) = 0. (60) 

Eq. (f5"9"f has solutions: 

• for a < 1: p 2 — 0, pi — 1 — a; 

• for a > 1: pi = 0, p 2 = 2=1. 

Eq. ([6H|l has solution only for p 2 — and pi = 1 — 6. This means that if the densities of the particles are 
non zero there is no phase transition. 



3 As explained in j 10] when a = 1 one has to take a single contour integral around 1. 
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5 Conclusions 



In this paper we have studied an ASEP with two species of particles and different hopping rates. We 
have formulated the computation of the cumulants of the currents as an eigenvalue equation, and we 
have shown that this leads to an integrable (d la Yang-Baxter) transition matrix. This has allowed us to 
employ the formalism of the Algebraic Bethe Ansatz to solve the problem, by finding the Bethe Equations 
for an arbitrary number of particles of each specie. The analysis of the Bethe equations gives in principle 
all the cumulants of the currents. We found the exact formula for the velocity of the particles of type 
2, and computed its limit when the size of the system goes to infinity, keeping non-zero densities for the 
particles. We find this way that, when the densities are different from zero, there is no phase transition. 

Our work can be extended in different directions. First, we think it would be interesting to use the 
Bethe Equation we found, to compute the spectral gap as a function of the hopping parameters a and 
/?. We have briefly discussed in the Appendix the extension of the problem to a larger number of species 
and different hopping rates, it would be also nice to work out the average velocity of particles of a given 
type as functions of the hopping parameters. Another interesting possibility is to consider the problem 
on a lattice with open ends and letting particles flow in and out of the system. We plan to come back to 
these issues soon. 
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A Algebraic Bethe Ansatz for the ASEP with one specie 

Let us consider the ASEP consisting of m particles and r holes on a ring of size N — p + r. Each particle 
can jump into a neighbor site only if the site is empty. The probability for jumping forward is q ■ dt, while 
the probability for jumping backward is p ■ dt. Following [9] we consider the total distance covered by all 
the particles in a time t, denoted by Y t . In order to determine the behavior of Y t we look at the joint 
probability P t (C,Y) of being at time t in a configuration C and having all the particles covered a total 
distance Y t = Y. The generating function 

oo 

^(C) = ^e^ r P t (C,y), 

Y=0 

which satisfies the following evolution equation 

^-F t (C) = Y,i M o(C,C) + e^M.iCC') + e'^M^(C,C')}F t (C). (61) 
at c , 

where M\{C,C')dt is the transition probability for going from the configuration C to the configuration C 
and moving a particle forward of one step, while M_i(C, C')dt correspond to a particle moving backward 
of one step and Mq is the diagonal part. The large time behavior of (e l ' 12Yt ) is determined by the largest 
eigenvalue \{u X2 ) of the matrix transition matrix M (C,C) + e Ul2 M X {C , C) + e - " 13 M_i(C,C'). 

The transition matrix M (C,C) + e Vl2 M 1 {C,C) + e~ Vl2 M-i(C,C) can be diagonalized by means of 
the Algebraic Bethe Ansatz [19, 20 . In our case we are led to consider 

Ra, b (x, y) = 1 + X(x, y)E (62) 
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where 



\(x,y) 



& — y 

e 2 



e 2 



p e 2 — q e 2 

The matrix R a ,b(x, y) acts on F a <g> H ->■ V a ® V & , where V = C 2 , and in the basis |00), |01), 1 10), |11) the 
matrix E reads 

/ \ 



E 







-q pe 







(63) 



qe" 12 -p 
\ J 

We define the matrices L a _ ai (x,yi) = Pa.atRa^iiXjyi), where P a ^ is the permutation operator. We 
introduce also a matrix f2 which acts only on the auxiliary space, in a diagonal way 



wi(j7) o 

u) 2 (y) 



and its entries can depend on some auxiliary parameters y. The reason for considering such a generaliza- 
tion of the problem we started from, which in facts correspond to fl equal to the identity, comes from the 
ASEP with two species, as seen in the text. For a system of size N the monodromy matrix is constructed 
by means of L a , ai (x, yi) 

^a®j^{x,y) = L a , afi {x,yf,) . . . L a>ai (x, yi). 



The transfer matrix is given by 



T (x,y,y) = tr„ ^l{y).%®^{y) 



The fact that [O(y) <S> £l(y), R(x, x')] = 0, combined with the Yang-Baxter equation, implies that the 
transfer matrices with different parameters x and x' commute among themselves. The transition matrix 
of the ASEP is obtained as the logarithmic derivative of the transfer matrix in zero (at SI = Id and 

Vi = o). 

Let us write the monodromy matrix in the auxiliary space as 



ZT{x,y) = 



Mx,y) B(x,y) 
C(x,y) D(x,y) 



(64) 



the Algebraic Bethe Ansatz proceeds by constructing an eigenvector acting with B(Q,y) on a reference 



state. Our reference state is |1) 




corresponding to the completely full system, 



which is an eigenvector of the transfer matrix. Indeed we notice that A(x,y)\l) = |1), D(x,y)\l) = 
n£=i (^F 1 ) |1>. 0(x,§)\l) = 0. Hence the eigenvalue of |1) is Wi(tf)A(x, §) +uj 2 (ff)D(x, §) = ^(y) + 

u>2(y) IlfcLi (p^( x i Vk)) ■ We search now for eigenvectors of the form 

|(Cl,...,Cr)) = B(Cl)-..-B(Cr)|l> (65) 

and use the Yang-Baxter algebra, satisfied by the operators A(x),B(x),C(x),D(x) as a consequence of 



Notice that ys and y can and will be in general different quantities. 
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the Yang-Baxter equation 



[A(x),A(y)] = 0, (66) 
A(x)B(z) = ^B(^)-^«B(^ W (67) 



B(z)B{x) = B(x)B(z) 

e Ul2 p" ia fl — n\(r z)) 

D(x)B(z) = ^—B(z)D(x) 1 q Y " B(x)D(z). (69) 

p\(x,z) pA(x,z) 

The requirement |(£i, ■ ■ • , Cr)) to be an eigenvector can be expressed in terms of a Bethe equation 

^n(|g)nww^ m 

which fixes the values of £. The eigenvalue of the transfer matrix is given by 



k=l 



We consider now the limit p — > 0, q — 1. In order to do that without getting a singular limit we have 
to change the spectral parameters into Q — > Q — l ° s ^ pS> . Then the Bethe equations turn into the form 

11 ^ p e (0-0)/2 _ g 6(0-0/2 J 11 ^ e (Cf-Sfc)/2 _ qe (vh-Ci)/2 J e [ U > 

and for p — > and q — 1 we get 



Defining Zj = e ^ we arrive at 



Sin(-|)n(^)-^ 

and the eigenvalue can be simply written as 

r 

A(x) = e^ 1 (y)rj(l-e^ l ) (75) 

i=l 

B Yang-Baxter equation for multi-species ASEP with different 
rates 

In this appendix we want to discuss to what extent the Baxterized form of the R— matrix ([7|) can be 
generalized, in order to describe a process with a number of species greater than 3, and a hierarchical 
structure. Labelling the species with numbers from 1 to n, the hierarchy means that a particle of kind 
i can hop to the right of a particle of kind j only if i < j. It is well known that if all these elementary 
processes happen with the same rate, then the R— matrix is simply given by the Baxterization of the 
Hecke algebra [18] . What we want to consider here is the case when the hoppings among the particles 
depend on the species involved in the hopping 
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ij ~ * ji with rate a,y if i < j. 

a/3,7<r 



Writing the matrix describing the hopping ij — > j'z as E^f) — S a idpj(e Uij 5 aa Sf} 7 — 5 a7 Sp a ), we would 



like to find an P— matrix of the form 

R{x,y) = l + Y,9ii(x,y)E {ij \ (76) 

which satisfies the Yang-Baxter equation 

Ra.biv, z)Rb, c {x, z)R a>b (x, y) = R b , c {x, y)R a , b {x, z)R b . c (y, z), (77) 

and such that 

dX x=y=0 f— ' 

im 

Let us define for i < j < k the projectors Pij and Pi,j,k, which act on C n (g>C™ (g>C n . Py projects on the 
states occupied only by particles of type i and j, while Pij,k projects on states occupied only by particles 
i, j and k. If we intertwine eq. (|77|) with Pjj we see that we reduce to the problem with two types of 
particles (or equivalently one type of particles and the empty sites), and we recover easily that g% must 
be of the form 

If we intertwine ea.(|77p with Pj,j,fc we recover the problem with three types of particles treated in the 
main body of this paper, and the Yang-Baxter equation implies that 

f i;k (x) = hjtx) + % frl{x) = fr£(x) + 6*;*. (80) 

This means that all the functions fij(x) are determined in terms of a reference one, which we chose to 
be fi } 2(x), and of the parameters b\\. and b l - k . Actually the relations in (f80|) put also constraints on the 
bs. Indeed it is easy to see that we must have 

fc-i 3-1 

°i.k ~ / ;i.l+V °i.k ~ / ,l+l.k- 

l=j l=i 

Moreover if i > 1 one can get fij+±(x) starting from/i.j(x) in two different ways 

fi,j — * fi,j+l or fi,3 ~^ fi—l,j ~ * + l ~ * fi,j+l- 

The previous relation fixes 

_ , h i,j _ °i-l,3+l 

°i.3'+l ~ 1 _ h i-h3 h i-h3 anQ °<.3+l - , _ h i-l,j ■ 

In conclusion we can chose as free parameters and 6;Zi';j which in a problem with n species 

are in number of 2(n — 1). Hence among the n(n — l)/2 rates ay, only 2(n — 1) are independent, given 
the form of the R— matrix l[7^|) . This of course does not rule out the possibility that the problem with 
generic rates is integrable, but one should look for a more general R— matrix to prove it. 
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